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Abstract 

We look at the transition to the semiclassical behaviour and the 
decoherence process for the inhomogeneous perturbations in the in- 
flationary universe. Two different decoherence mechanisms appear: 
one dynamical, accompanied with a negligible, if at all, entropy gain, 
and the other, effectively irreversible dephasing, due to a rapid vari- 
ation in time of the off-diagonal density matrix elements in the post- 
inflationary epoch. We thus settle the discrepancies in the entropy 
content of perturbations evaluated by different authors. 



1 Introduction 

The understanding of the origin of the large-scale structure in the universe 
remains one of the central problems in modern cosmology. It is believed 
that structures arise because the universe in the early epoch was not exactly 
homogeneous and isotropic, but must have contained some density irregular- 
ities. 

These density irregularities or fluctuations may be classified into two ma- 
jor classes. The first one of these, the primordial fluctuations, corresponds 
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to those present in the universe ab initio and which may be thought of as 
classical irregularities in the initial structure of the universe, as for example 
the Mixmaster universe [I], the similarity solutions [EJ, gravitational solitary 
waves ||, etc. 

A more exciting possibility, and the one relevant to this paper, are the dy- 
namically induced fluctuations, those arising inevitably in exact FRW models 
due, for example, to graviton creation in a time varying background gravita- 
tional field M. These fluctuations are "seen" as inhomogeneous fluctuations 
of the geometry. 

Since the primordial fluctuations are classical in their nature no prob- 
lem arises in interpreting the density or temperature variations. The prob- 
lem with the dynamically induced quantum fluctuations is more involved, 
however, and touches one of the most fundamental puzzles of physics, the 
decoherence process in the process of quantum to classical transition. 

In a recent series of interesting papers Lesgourgues, Polarski and 

Starobinsky (LPS in what follows) analyze the evolution of inhomogeneous 
perturbations generated from the vacuum state during the initial accelerated 
expansion of the universe. These authors find that the quantum fluctuations 
become classical with stochastic Gaussian amplitudes, and that the decaying 
mode of these fluctuations becomes exponentially small towards the end of 
the inflationary period. Consequently, LPS suggest that the exponentially 
small decaying mode may be discarded for all practical purposes in this stage 
of the evolution, and the initially quantum perturbations may be considered 
classical. It seems then that there should be no (or minimal) information 
loss, and, apparently no (or minimal) entropy generation in this quantum to 
classical transition. 

It is usually believed that this transition occurs when some class of coarse 
graining is enforced on the system. LPS question the validity of different 
coarse graining schemes in the following sense: if the dynamics of the system 
leads in a unitary way to a semiclassical behaviour is it necessary to per- 
form a coarse graining? Last, but not least, stands the fact that the results 
obtained by LPS on the entropy content of the perturbations are in strong 
disagreement with respect to the results obtained by other authors (|| and 
references therein). 

The study of the entropy of quantum fields in the cosmological context 
dates from the early days of the developments in quantum field theory in 
curved spacetime (see || for a review up to the 80's), with further elaboration 
by Hu, Kandrup, and others |10| . Currently the language used to discuss the 
subject is that of squeezed states, introduced in this context by Grishchuk 
and Sidorov [II], and borrowed from the quantum optics community [12 



The squeeze formalism was also recently used by two of us to study the 
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problem of quantum tunneling under the influence of a time varying force, 
and the speculation was put forward that the initial entropy gained by the 
universe in the "creation from nothing" picture may be simply evaluated by 



using the squeeze parameter [13 



The use of the squeezing formalism should not obscure the fact that, 
in the standard view, the entropy growth is always due to the particular 
coarse graining one chooses; yet, interestingly enough, the whole variety of 
coarse graining procedures (averaging over a period of the squeeze angle 



13 1 , integrating over rotation angles [fL4|, neglecting information about the 



subfluctuant variable ]T5|], setting off-diagonal elements of the density matrix 



to zero [0) leads in the large squeezing limit to the same expression for the 
entropy generated per mode, S& ~ 2r^. LPS M, however, predict almost no 
entropy generation for similar systems. 

One of the main purposes of this paper is to show, within the example 
considered by LPS, that as long as the Wigner function is used to represent 
the evolution of the system, as they do, the dynamical limit (certain pa- 
rameter £ — > oo in our parametrization) and the semiclassical limit, % — > 0, 
are identical, and probably for most practical needs the system they con- 
sider may be labeled as "dynamically semiclassical". However, the delicate 
question of the entropy content is related to the limiting behaviour of some 
elements of the density matrix associated with the system. We show that for 
the density matrix the £ —>■ oo and H — > limits are different, and that the 
"true" decoherence is achieved only in h — > limit, due to effective dephasing 
or rapid oscillations of the off-diagonal elements of the density matrix. This 
we do in the following Section. 

In Section 3 we consider a model universe which starts inflating and then 
passes to radiation and matter dominated epochs. Evaluating the expression 
for the entropy based on the density matrix calculations we find that, due to 
an effective dephasing, the entropy generated per mode is given by the usual 
expression 6\. « 2r k . 

Section 4 is dedicated to some final remarks. 



2 Wigner Function and Density Matrix 

The study of the cosmological perturbations generated after the amplification 
of the vacuum fluctuations in the early stages of the universe can be reduced 



to the analysis of the evolution of a scalar field in a FRW background |T7 . 
If at some initial conformal time rjo the field is in the vacuum state, the 
dynamical evolution under the influence of external gravitational field drives 
the scalar field into excited energetic states with opposite momenta. 
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One way to analyze the dynamical evolution of the cosmological per- 
turbations and their quantum-to-classical transition is through the Wigner 
function formalism. The Wigner function for a one-dimensional quantum 
mechanical system is defined by ITT 



W?,P)~ J^dx (q-x/2\p\q + x/2) (1) 

where p is the density matrix of a pure or mixed quantum state. 

The Wigner function corresponding to the evolution of the vacuum state 
of a scalar field in a FRW background was evaluated by Polarski and Starobin- 
sky P, and is given by 



W(k , _ k) = _^ exp |_Mp!} exp {-\f W k) - p„(k)r \ , (2) 



where y(k) is the Fourier transform of the rescaled scalar field (y = a<ft) 
and p(k) is its canonical conjugate momentum. F(k) and \fk\ are two time- 
dependent parameters related to the variances of the field and the momen- 
tum. They are related to the squeeze parameter and the squeeze angle 
(Pk, the characteristic parameters of the Schrodinger picture, through || 

\fk\ 2 = ^(cosh2r fc + cos2v9 fc sinh2r fe ) F(k) = -sin2<£> fe sinh2r fc . (3) 

If we decompose the field y(k) and the canonical conjugate momentum 
p(k) into their real and imaginary parts, the two dimensional Wigner func- 
tion may be further expressed as a product of two identical one- dimensional 
Wigner functions 

W(k, -k) = Wi(k, -k)W 2 (k, -k), (4) 

where W\ and W2 are related to the real and imaginary part of the field and 
momentum. 

Both W\ and W2 may be thought of as a particular case of the general 
parametrization of the Gaussian Wigner function (see for example Cooper et 
al. @): 

„2 ot2 / \ 2 




(5) 

From now on we will use the expression fl5|) for convenience. 

Here, comparing with Eq. (0), the parameters £ = V2\fk\ and p = 
F(/c)/(v / 2|/ fe |) are related to the variances of the perturbations of the field 
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and initially correspond to a coherent state. Note, as well, that this Wigner 
function is identical to the one obtained in the problem of an upside-down 
harmonic oscillator [ |2U[| . 

The Wigner function, in fact, is only another (particularly adequate) way 
of writing all the elements of the density matrix. The relation is given by 
the following Fourier transform 

(x'\p\x) = [ +0 °dpe^'-^ h f w (?-^,p). (6) 

J -co Z 

Using the expression (M), we obtain for the density matrix the following 



expression (cf. Cooper et al. fl9|j) 



(APMM = j^yT 2 exp {-^' 2 + * 2 ) + |^' 2 - * 2 )} • (7) 

It is clear that this density matrix corresponds to a pure quantum state, for 
it satisfies p 2 = p. 

Let us analyze the behaviour of the Wigner function in two different 
limits. One, which we will call, following ||, the dynamical limit and the 
other the semiclassical one to be defined below. 

Since we are interested in analyzing the dynamical evolution of cosmo- 
logical perturbations, we consider a limit corresponding to £ — > oo along 
with p/l; being kept constant. This limit is related to the large squeezing 
limit characterized by the behaviour of the parameters F(k) and \fk\ which 
become unbounded. Thus, the dynamical limit obtained from the Wigner 
function (|5|) is 

/ H V(x,p) = ^=e-^ 2 ( 5(p-fx). (8) 

On the other hand, we can express the semiclassical limit in the usual way 
by analyzing the small h expansion in the Wigner function ([5]). It follows 
that 

ffi{x,p) = &{x,p)- (9) 

Both limits represent a classical probability distribution in the phase space 
with x obeying a Gaussian law, whereas p is fixed by the value of x at any 
instant. 

Since both limits give the same expression for the Wigner function, it 
seems that the dynamical limit and the semiclassical limit are equivalent. 
Thus, for any practical purpose (physical measurement of different magni- 
tudes) one expects to obtain the same result. 

Yet, these two limits are not completely equivalent. If one tries, say, to 
reconstruct the elements of the density matrix using the equation (|6|), one 
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then immediately runs into trouble with the expression for f^(x,p), since it 
is valid only for small h, and, therefore, one can not use the Fourier transform 
involving e ip( - x ~ x )/ n ^ given that it is not perturbative in h. On the other hand, 
this Fourier transform is perfectly defined for the expression f^{x,p) which 
was obtained for large £ instead. 

Let us reinforce the statement by looking at the behaviour of the system 
following these two limits but starting directly from the density matrix (0). 

The dynamical density matrix pd, i.e., the dynamical limit of the density 
matrix as £ — > oo, is precisely the one given by the expression ©, in the 
sense that pd = p, even for very large £. The unitarity of the evolution of the 
system is obvious in this limit and pd satisfies Trp^ = 1. 

The leading term of the semiclassical expansion (U — > 0) of the density 
matrix corresponds to 



\3' \Psc\Xj 



2£h 



6(x)8(x'). (10) 




We can see that the semiclassical limit and the dynamical limit for the density 
matrix are quite different: pd ^ p S c, with the difference between them being 
in subleading terms in fi. 

One may still define yet a different limiting density matrix obtained from 
the dynamical limit of the Wigner function since it is quite "legit- 

imate" as explained above. Performing the Fourier transform (|6|) we get 

/ l\~ | \ 1 \ W + x) 2 ip /2 

<?\p*\*) = (^a^ exp { — s^ + 2^ (x 

= (x'\p d \x)exp^- ( X ^ | . (11) 

We see that the non-diagonal elements of the density matrix obtained by 
Fourier transforming the dynamical limit of the Wigner function /^(x,p) dif- 
fer from those of the pure dynamical density matrix pd in subleading terms. 
Furthermore, the density matrix ( |TTD does not satisfy p 2 d = pd and gives 
Trp^ = oo. This entails that pd neither represents a pure state, nor corre- 
sponds to a quantum density matrix. 

We now look at the entropy. The quantum evolution of the system cor- 
responding to pd is unitary (it is governed by a quadratic time- dependent 
Hamiltonian) . Therefore, there is no entropy change associated with the 
quantum evolution of the system in time, and since the initial state is pure, 
it stays pure forever giving zero von Neumann entropy 

S d = -Trpd\np d = 0. (12) 
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Let us now turn to the density matrix pd obtained from the dynamical 
limit of the Wigner function. Since pd differs from pd in subleading off- 
diagonal terms, the "entropy" that were to be associated with it should 
tell us how much information would be lost by discarding those subleading 
terms, or, in other words, how strongly "mixed" the state given by the density 
matrix pd is, as compared to the pure state given by pd- 

However, the von Neumann entropy — Trp^ In pd is not well defined here. 
The fact that Trp^ = oo and that Trp^ = 1 implies that the matrix has 
negative eigenvalues. It follows that different branches of the logarithm would 
give different imaginary contributions to the von Neumann entropy, if defined 
in the usual manner. We tried the following regularization: 

Sd = -Tr-Pdlnpl (13) 

Intuitively, one might think that the quantity just defined is the entropy 
gain due to a dynamical decoherence. If so, it must be consistent with the 
calculations of LPS. However, one may easily show that this entropy has the 
wrong sign and is clearly pathological. The telling point here is that the 
state described by pd is ill-defined, and the calculations produce pathological 
results. This shows that one must be extremely careful when such a delicate 
quantity as the entropy is considered. 

The dynamical limit of the density matrix pd given by expression (|7|) is 
the one which correctly describes the quantum state of the system. In fact, 
the system remains quantum forever. To obtain a classical behaviour one 
should in principle discard the quantum interference effects. These effects 
are associated with those terms of the density matrix which oscillate rapidly 
due to the smallness of h. Thus, by averaging the quantum pd matrix over 
classical observation times (which is feasible in the cosmological context), 
the off-diagonal terms in the position representation are washed away. When 
the same procedure is applied to pd, the density matrix obtained from the 
dynamical limit of the Wigner function (||) , we get exactly the same averaged 
density matrix, since both density matrices have identical rapidly oscillating 
phases given by the expression exp \^^(x' 2 — x 2 )|. 

Performing the average we have just mentioned, one readily finds an ef- 
fective coarse grained entropy 

S^ln£, (14) 

or as expressed in terms of the squeezing parameter r, (£ ~ cosh2r), the 
entropy per mode is 

S^2r, (15) 
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in the large r limit. 

This result is the same as the one we obtain by evaluating the Boltzmann 
entropy from the probability in the phase space described by the Wigner func- 
tion fyy(x,p) after averaging over the momentum f^(x) = Jf^, f^-(x,p)dp. 
The above coarse graining scheme is thus equivalent to tracing out the non- 
diagonal elements of the density matrix and keeping the diagonal elements 
alone. 

In order to clarify the issues raised by the "bad" behaviour of the elements 
of the matrix density in the position representation, we shall work in another 
representation which will prove to be more convenient. 



3 Dynamical Evolution of the Density Matrix 

As previously pointed out, we are interested in the way quantum interference 
effects can be suppressed following the dynamical evolution of the system. 
The behaviour of the system at large times, i.e., when the inflationary epoch 
is followed by radiation and matter dominated periods, is similar to that of 
a parametric oscillator whose time-dependent potential is slowly turned off. 
Thus, it seems that the most adequate basis to study the evolution would be 
that of the harmonic oscillator. 

In the Schrodinger representation the dynamics of the state is given by 
the time evolution operator |0, r])s = S\0, 770) - The time operator S is the 
usual two-mode squeeze operator and the time evolution of the density matrix 
expressed in the harmonic oscillator basis is given by |12|, |19| 



(n = 2l\p\m = 21') = (n\S(r, (p)\0)(0\S*(r, (p)\m) 



, ((2/)!(2Q!)V 2 (tanhr)^ itp{l _ V) 

2 l + l 'l\l'\ coshr ' 1 ' 

for each mode k. For convenience we have dropped the sub-index k, and r 
and cp are the usual squeeze parameter and squeeze angle, respectively. The 
only non vanishing diagonal elements are 

/9/U9A (2/-l)H (tanhrf 

<2J|p|2i> = (17) 

Note that only the non-diagonal elements of the density matrix depend on 
the squeeze angle ip. 

The time dependence of the density matrix elements in this basis is given 
by the evolution of the squeezing parameters which can be obtained through 
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the interrelation between the Heisenberg and Schrodinger pictures [T^, |2l 
The equations of motion for r, ip and 9 are H 12), ^1 



. a 

r = — cos 2^9, 
a 

= -fc - — coth2r sin2y?, (18) 
a 

a' 

0' = +k H tanhr sin2y9. 

a 

The squeeze parameter r grows indefinitely towards the end of the inflation- 
ary period while the squeeze angle <p approaches a constant value. Thus, the 
non-diagonal elements of pd do not oscillate, and no dephasing takes place, 
as was correctly pointed out by LPS. 

However, let us now push the evolution of the system forward, pass the 
inflationary period by matching it to a radiation dominated one followed by 
a matter dominated expansion as, for example, was considered by Grishchuk 
and Sidorov jTJ]. It follows then that the behaviour of the squeeze parameter 
r and the squeeze angle <p are qualitatively different from the one we see 
during the end of the inflationary period. 



To be specific, we use Grishchuk's and Sidorov's |TI]] simple model of the 
universe containing the three stages of expansion just described above 

^(v-ZVi) (Vi<V<Vz), (19) 



-<x> < rj < rji < 0), 



(r) - 4r)t + 7] 2 ) 2 (t] 2 < T] < oo) 



where at, a r and a m represent the scale factor during the inflationary era, 
radiation dominated period and matter dominated period respectively, and 
rji and i]2 are times where the transition between two consecutive periods 
takes place. 



We have solved numerically equations (fL8|) under conditions (|T^) and 
found the following behaviour of the squeeze parameters: the squeeze param- 
eter r approaches its maximum value at the end of the inflationary period, 
oscillates around this large maximum value during the radiation dominated 
period and finally sets to a constant in the matter dominated epoch. This 
behaviour can be seen in the Figure 1, where the parameter r is represented 
as a function of conformal time rj. 

On the other hand, as is apparent from Eq. (fL6|), the squeeze angle <p 
shows itself as an angular phase e 21 ^ 1 - 1 ') in the non-diagonal elements of 
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Figure 1: Dependence in time of the squeeze parameter r. The following 
values for the parameters appearing in the dynamical equations are chosen: 
7] = —20, 771 = — 10~ 5 , 7] 2 = 20 and the wave number has been set to k — 1. 



the density matrix. This angular parameter (p grows linearly with time as 
tp ~ ifQ — kri leading to oscillating behaviour of the off-diagonal density matrix 
elements. Thus, the dynamical evolution of the model suggests a natural way 
to average these elements to zero. It is clear then, that sooner or later an 
effective dephasing process will take place. A typical behaviour of the phases 
of the off-diagonal density matrix elements is presented in Figure 2. 



Cos cp Sin 9 
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Fig. 2a Fig. 2b 

Figure 2: The figure 2a represents the real part of the phase of the off-diagonal 
density matrix elements. Imaginary part of the phase is represented on figure 
2b. 

As a consequence of the above dephasing process one may proceed in cal- 
culating the entropy in the usual way by dropping the non-diagonal elements 
of the density matrix. 

00 

S = -]T(2Z|p|2/) (2/|lnp|2/) 
1=0 
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~ (21 - 1)!! (tanhr) 2 ' , 
2^ — — : In 



1=0 



2 l l\ 



coshr 



(2/ - 1)!! (tanhr) 



21 



2 l l\ 



coshr 



(20) 



For convenience, the latter expression can be separated into two pieces 
s g (21 - 1)!! (tanhr) 2 ' ^ [(tanhr^ 2r 



1=0 



2H\ 



coshr 



and 



S- 2 



~ (2/- 1)!! (tanhr) 2 ' ^ 



2H\ 



coshr 



coshr 



(2Z-1)!! 
2H\ 



where the first one readily gives 

Si = ln(coshr) — sinh 2 (r) ln(tanhr). 



(21) 



(22) 



(23) 



To evaluate the second term we expand expression ( ]22|) in terms of gamma 
functions and get for r ^> 1 



So 



E 



r(2Q 



coshr ^ 2 2 '-!/(r(/)) 2 



(tanhr) In 



F(2l) 



2 2l - i i(v(i)y 



(24) 



With the aid of the Stirling's formula In 

1 



r(2Q 



2 2i - 1 z(r(0) 2 



-~ ln(/7r), we obtain 



S, 



2 cosh r 



« (tanhrf , 
— w^— ln ( /7r )- 



(25) 



Using the McLaurin asymptotic method [22| for large r we get 



So 



1 



1 



2y/2 coshr a/— In tanh r 
+ Of— In tanhr], 



7 + log(- 



-8 In tanhr, 



c 



71 



(26) 



where C is a constant. 

Finally, the total coarse grained entropy is given by 



S 



2\[2 coshr 
and in the limit r — > oo we obtain 



In cosh r — sinh r In tanh r + C — 

1 (— In tanhr) _1 / 2 / —8 In tanhr, 

- 7 + log(- 



7T 



(27) 



S«2r. 



(2f 
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4 Conclusions 



Before closing we feel that some final remarks are in order. 

The first lesson learnt from our analysis is that one must be careful 
with the equivalent expressions obtained within different limiting procedures. 
From the example with the Wigner function we have seen that although the 
two expressions coincide (the one obtained in the dynamical limit, £ — > oo, 
of the Wigner function and the other in the semiclassical one, h — > 0), one 
cannot use the last one to reconstruct the elements of the density matrix. 
The reason is simple: the reconstruction of the density matrix involves again 
the parameter %. For any other calculation, these two limits may be perfectly 
equivalent, but if one is interested in the entropy of the system, the procedure 
fails! One may easily see this also by evaluating the two limits starting di- 
rectly from the expression for the density matrix elements (|7|). The difference 
between the limits is evident from expression fllCf ). 

Going somewhat further, we were interested in seeing whether, by com- 
paring the elements of the density matrix of the pure state taken in the late 
time limit and those "legitimately" obtained from the dynamical limit of the 
Wigner function, we can get a measure of the information loss due to dy- 
namical decoherence. Intuitively, this measure should have been consistent 
with the small gain in the entropy due to the phase-space volume increase 
as described by LPS @. Unfortunately the entropy expression calculated in 
this way is ill-defined. 

To evaluate properly the entropy generation for each mode of the pertur- 
bation we study the time evolution of the density matrix in the harmonic 
oscillator basis. Physically, this seems to be the correct basis since at late 
times, during the post-inflationary behaviour, the time dependent potential 
governing the dynamics of the modes switches off adiabatically. We explicitly 
see the dephasing behaviour due to the rapid oscillations of the off-diagonal 
elements of the density matrix, and calculate the entropy growth recovering 
the result Sk ~ 2rfc. It is important to note, however, that each mode has its 
proper time scale for dephasing depending on its wavenumber k, yet (almost) 
all modes, as mentioned above, sooner or later contribute to the process. 

Decoherence is normally understood as the vanishing of the off-diagonal 
density matrix elements in a particular basis and is accompanied by the 
entropy growth due to the loss of information encoded in those elements. We 
have evaluated this entropy and found it consistent with previous results. 

On the other hand, we have found that the entropy associated with the 
dynamical decoherence process evaluated on the basis of the dynamical limit 
of the Wigner function is ill-defined and probably carries no useful informa- 
tion about the statistical properties of the system. 
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